Time domain analysis of dynamical switching in a Josephson junction 
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^SJ ' We have studied the switching behaviour of a small capacitance Josephson junction both in experi- 

ment, and by numerical simulation of a model circuit. The switching is a complex process involving 
the transition between two dynamical states of the non-linear circuit, arising from a frequency de- 
pendent damping of the Josephson junction. We show how a specific type of bias pulse-and-hold, 
can result in a fast detection of switching, even when the measurement bandwidth of the junction 
voltage is severely limited, and/or the level of the switching current is rather low. 
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r-H ■ I. INTRODUCTION 

D \ Josephson junction physics is presently experiencing a renaissance, where many new and exciting experiments on 
G . small capacitance Josephson junction circuits are demonstrating quantum dynamics with measurements in the time 
domain iSi^iii^ In this field, an engineering approach to the design of quantum circuits is presently being applied to 
C^ , the long term goal of realizing a solid-state, scalable technology for implementation of a quantum bit processor. The 
intrest in quantum computation is fueled by a theoretical dream of massively parallel computation with quantum 
two-level systems, or qubits (for a nice introduction, see the book by Nielsen and Chuang®). A central question for 
■ O ■ these Josephson junction experiments relates to the optimal design for qubit readout, or detection of the quantum 
state of the circuit. Quantum state readout has been succesfuUy performed by switching current measurements. 

In this paper we use numerical simulation to analyze a particular type of readout -the switching of a Josephson 
junction subject to frequency dependent damping. The system we model is motivated by experiments, and comparisons 
are made with measured data. The switching process is a complex transition between two dynamical states of the 
circuit, and our numerical simulations allow us to study the speed and resolution of the switching process, so that we 
can investigate the use of this switching process as a detector for quantum state readout of a Cooper Pair Transistor 
(CPT). 

l/~j The classical description of the Josephson junction reduces to an analysis of the non- linear dynamics of the Josephson 

\^ I phase variable, which is conveniently discussed in terms of a ficticious "phase particle" moving in a washboard or 
^D ' cos((/)) Josephson potential.^ The time evolution of the phase variable depends entirely upon the particular circuit in 
\l . which the junction is embedded. The simple and most tractable model is the so-called resistively and capacitively 
^^ ■ shunted junction model (RCSJ) of Stewart^ and McCumber^, where an ideal Josephson tunneling element is connected 
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in parallel with an ideal capacitor (the tunnel junction capacitance), an ideal frequency independent impedance (a 

resistor), and driven by an external current source. More difficult to describe, and more relevant to experiments 

H \ with small capacitance Josephson junctions, is the phase dynamics when the parallel impedance becomes frequency 

i-rt dependent, and when fluctuations due to the finite temperature of the dissipative elements are included. With small 

(~| ' capacitance Josephson junctions we are often in a situation where the frequency dependent shunting impedance gives 

O 1 overdamped phase dynamics at high frequency, and underdamped phase dynamics at low frequency. The generic 

J-^ ' features of such a model have been studied extensively with computer simulations by e.g. Kautz and Martinis.^" In 

^ , this model the system can be in two distinct dynamical states for some given bias currents. In the "phase diffusion" 

! ""j • state, a small finite voltage appears over the junction when the noise term causes a diffusive motion of the phase 

rS particle in a tilted washboard potential. In the "free running" state, a much larger voltage is sustained across the 

jrt ' junction. A switching of the junction, or transition between these two dynamical states, occurs when a critical velocity 

■ ■ ■ ' of the phase particle is reached, so that it can overcome a "dissipation barrier" resulting from the frequency dependent 

dampingJiiSiiiiiS. In the present work we study the probability for this dynamical switching, as a function of time 

and amplitude, for a particular type of bias pulse. 

The bias pulse, which is shown in Fig.^a), consists of a short switch pulse of amplitude Ip and duration Tp, followed 
by a long hold level with amplitude Ihoid and duration Thoid- The idea behind this pulse is to quickly accelerate the 
phase particle above the critical velocity for dynamical switching, in order to force the switching process to happen 
as quickly as possible. The hold level is used to maintain the free running state with high voltage, long enough so 
that the voltage can be measured with an amplifier at room temperature. The hold level must be set high enough 
so that the phase particle will not be retrapped in the phase diffusion state, and yet low enough not to induce a 
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FIG. 1: (a) The bias pulse and hold, and (b) the frequency dependence of the shunt impedance Re[Z] (dashed line) and of the 
quality factor Q (solid line). 



late switching event. This type of bias pulse exploits the latching nature of the circuit to realize a sample- and- hold 
measurement strategy, in order to overcome measurement bandwidth limitations resulting from filters, lossy cables, 
and high gain amplifiers. The latch gives us the added advantage of having a binary detection of switching, with two 
distinct voltage outputs for switch or no-switch, from one bias pulse. 

The frequency dependent damping used in our numerical simulations is shown in Fig. ^b) , where we have plotted 
the real impedance Re[Z] of the bias circuit (dashed line), and the Q factor of the junction (solid line) as function 
of frequency (see sec. HTl and IIIII for junction and model parameters). The model used for the shunting impedence is 
a simple RC shunt (see sec. IIII|I which reproduces the qualitative feature of the actual experimental situation-that 
of overdamped phase dynamics at high frequency {Q <^ 1), crossing over to underdamped phase dynamics at low 
frequency {Q :$> 1). The frequency scale associated with the switch pulse (l/rp) can in principle be adjusted so that 
the junction is overdamped, and no hysteresis or latching behaviour of the circuit would result for the switch pulse 
alone. However, the frequency scale of the hold pulse (l/r^oid) is such that underdamped dynamics is realized, and a 
hysteresis or latching of the junction voltage is possible. 



II. THE EXPERIMENT 



The sample used in the experiments is actually a superconducting quantum interference device (SQUID), which acts 
as an effective single junction and allows us to continuously tune the critical current, /q, from its maximum value to 
zero by applying a magnetic field perpendicular to the loop. Here, however, we have only simulated data taken at 
a single value of Iq. The SQUID was fabricated by two angle evaporation of Al through a shadow mask which was 
defined by e-beam lithography. A micrograph of the sample is shown in the inset of Fig.|2Ia) together with a schematic 
of the experimental set up. An arbitrary waveform generator (AWG) biases the sample through an attenuator, a bias 
resistor via twisted pair cables of konstantan wire of length 2.1 m and DC resistance of 135 fi. The twisted pair acts 
as a lossy transmission line with a high frequency impedance of '^ 50 17. We have measured this value at frequencies 
up to 500 MHz. The losses in the line arc such that reflections or standing waves were only very weakly visible. At 
higher frequency the actual impedance seen by the junction depends critically on the geometry of the sample mount 
between the leads and the junction, which we have not measured or modeled here. No additional cold filters were 
added at low temperature. The low frequency impedance seen by the junction is given by the bias resistor, R^j ~ 2 kf7. 
The voltage across the sample is measured at the biasing end of the twisted pair cable through a gain 500 amplifier 
with bandwidth 100 kHz. The limited bandwidth of this amplifier does not allow us to observe the fast rise of the 
voltage when the junction switches to the free running state, the simulations however allow us to model the junction 
voltage in time. The bias resistor is chosen so that the voltage which builds up on the junction is kept well below the 
gap voltage V2A = 400 /ieV, in order to avoid quasiparticle tunneling in the junction and the associated dissipation. 
A digital counter registers a switching event when the output of the amplifier exceeds a trigger level. The switching 
probability is measured by counting the number of switching events and dividing by the number of applied bias pulses 
(typically 10^)^^ 

The capacitance of the SQUID junctions (parallel combination) is Cj ~ 7.2fF, as estimated from the measured 
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FIG. 2: (a) The experimental setup and sample shown in the inset, (b) The ideal version of the system consists of the Josephson 
junction (modelled by an ideal part and a capacitor Cj), connected to a current-source via some impedance Z{uj). (c) The 
impedance Z we use. Note that we have included the noise currents due to the resistors. 



junction area and the specific capacitance of 45 ^ / jim? . This gives a charging energy Ec — A.f? jlCj = 45 /ieV, or 
Ec/ks = 0.52 K. The normal state resistance of the junction is 1.18 kil giving a critical current /q = 265 nA and 
Josephson energy Ej = hIo/2e — 540 fieV, or Ej/kB = 6.3 K, or Ej /Ec ~ 12. This large ratio, together with the 
low tunneling impedance of the junction, means that we are justified in treating the phase as a classical variable. 
Measurements were made in a dilution refrigerator with base temperature of 25 mK, well below the superconducting 
transition temperature of Al, Tc ~ 1.2 K. The temperature is varied in the range 0.025 < T < 0.7 K as measured by 
a ruthenium oxide resistive thermometer in thermal equalibrium with the Cu sample mount. However, the effective 
temperature describing the noise in the theoretical model, is expected to be higher than the measured temperature 
as this noise is generated by dissipative elements located at higher temperatures. 

The AWG is limited to a rise-time Trise = 25 ns. We have measured that this rise time is effectively transmitted to 
the sample, although some of the switch pulse amplitude is lost for the shortest pulses due to dispersion in the twisted 
pairs. The switch pulse has magnitude Ip and length Tp (including 2Trise) followed by a hold level of magnitude 
Ihoid < Ip and duration Thoid- Every such bias pulse is followed by a wait time t^,, with 7 = 0, where the phase 
retraps and the junction returns to thermal equilibrium. 



III. THE MODEL 



A model circuit for the measurement set up, shown in Fig. |2Ib), consists of an ideal Josephson junction with the 
current and voltage given by the Josephson relations, in parallel with the ideal capacitance of the tunnel junction Cj. 
This ideal tunnel junction is biased by a linear circuit with impedance Z{uj) in parallel with an ideal current source. 
(The model could have equivalently been cast in terms of a series combination of an impdeance and voltage source). 
The model used for the impedance consists of a series combination of a resistor R2 = bOVt and capacitor C2 = 0.14 nF, 
in parallel with the bias resistor Rt, = 2 kH. (shown in Fig.l^c)). With these values we match the measured impedance 
of the twisted pair fairly well. This type of model, with high impedance for low frequencies and low impedance for 
high frequencies (Fig. ^b)), was first proposed by Ono et al^, and it is the perhaps simplest possible model of a 
frequency dependent environment that is known to correctly model the qualitative features of the switching process. 
The sample is at a low temperature (F = ksT/Ej in dimensionless form) and the resistors are at an unknown noise 
temperature, r„ — kgTn/Ej, which give rise to the noise currents /„b,2- 

Applying Kirchoff's rules and the Josephson relations to the circuit in FigE^b), one obtains a differential equation 
for (j){t). Defining the quantities t' = t/ts where tg — fi/2eIoZQ and Zq — Z{Q), to' ~ tsto, i = I /Iq (we will throughout 
this text use upper case / for currents in amps and lower case i for (dimensionless) currents in units of Iq) and 
Qo — Zoy^zeloCj/h, one can write the equation in dimensionless form. 
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where Ui{(t)) 



is the tilted washboard potential, z„ is the noise current due to Z with a correlation 



function obeying the fluctuation-dissipation theorem-'^^ (*n(^i)«n(^2)) = r„— /_ °° ReY{Lu') cosuj'{t[ — t'2)dw' , and the 
last term is the friction term in the form of a convolution integral between the voltage and the Fourier transform of 
the admittance, y{t') = J^^ Y{uj')e"^'^' diu' with Y = 1/Z. Note that if Z = i?, one recovers the RCSJ model. 

Using our specific model for Z, we can write Eq. ^ as three coupled first order differential equations-a form more 
suitable for numerical analysis. We find 
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where we have defined the following dimensionless parameters; 

V Vc RbCj 



^'c = -m-' P 



Rblo Rbh R2C2 



^^ l2eIoCj ^ 1 l2eIoCj 



h ' ^' l/Rb + l/R2 

where v is the dimensionless voltage over the junction and vc is the dimensionless voltage over C2. inj is the noise 
current due to resistor Rj (j = 6,2), with correlation function («nj(ii)inj(^2)) = 2r„^(5(i2 — i'l) where we for 

simplicity have assumed that Rb and R2 are at the same temperature. The generic features of the solutions to Eq. (O 
are known and have been investigated thoroughly by e.g. Kautz and Martinis^ 
The frequency dependent quality factor, plotted in Fig. ^b), is 

Q(-) = Q^^±Qkp-ML_ (3) 



where ujp = y^2eIo/fiCj is the plasma frequency (indicated in Fig.^b)). This expression is derived using the following 
definition; Q{uj) = ^ 2elQC j / hG^ {lj) where G = Re[l/Z]. 

We solve Eq. ^ using a 4'th order Runge-Kutta algorithm. The noise curents are inb,2 — Nab.2, where N G [—1,1] 

is a Gaussian distributed random number, ab.2 — a/ At^lT^ ^^ ^^^ standard deviation and At' is the dimensionless 
timestep in the numerical routine. This last factor is present since it is the mean value of inb,2 during a time At' that 
enters Eq. Q. 

We have simulated the switching process using Tp — 0.1, 1.0 and 10 fis. If the junction switches, there will be a 
finite voltage {V) > Vtrig over the junction (where {V) denotes the mean value over time and Vtrig is the trigger 
voltage) and the counter will click. The scale of (V) is determined by Rb and Ihoid- The finite voltage is latched for 
the duration of the hold pulse due to the hysteretic behaviour of the non-linear circuit. In the no-switch case, the 
phase particle either fluctuates thermally {{V) = 0) or diffuses down the washboard {{V) < Vtrig )-in both cases the 
counter will not click, assuming Vtrig is chosen high enough. 

IV. RESULTS AND DISCUSSION 

Two important quantities for the use of this switching process as a detector of the quantum state of a circuit, are the 
measurement time and the resolution. The measurement time can be determined by recording when the switching 
events occur, i.e. at what time the simulated voltage reaches the prescribed trigger level. Due to bandwidth limitations 
in the actual experiment, we can not measure this time. However, we can simulate many switching events to determine 
the probability, P{t)dt, that the switching event occurs between the time t and t + dt. A typical distribution P{t) is 
shown in Fig. |3fa). We define the measurement time to be the width of 98% of this distribution, neglecting the first 
and last percent. The resolution of the detector is determined from the probability, P(lp)dlp, that a switching event 
will occur with a switch pulse of amplitude Ip. This probability distribution can be both simulated and measured, and 
a typical curve is shown in Fig. IJ^b). The resolution of our detector A/, is the ability of the detector to discriminate 
between two different values of the switching current (which will differ depending on the quantum state of a circuit). 
We have (arbitrarily) chosen a discriminating power a — 0.8 (meaning that 20% of the events will be miscounted) to 
determine the resolution A/. 
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FIG. 3: (a) A probability distribution of switching as a function of time, with the measurement time defined as the width of 
98% of this distribution, (b) A graphical representation of the discriminating power a and the current resolution A/. 



Fig.^lshows a simulated Pit) for three different values of the hold level, ihoid- The magnitude of the noise was set 
to r„ — kBTn/Ej = 0.47 (corresponding to r„ = 3.0 K with Ej from secEJ and the voltage trigger to Vtrig — 35 ^V, 
which is safely above the phase diffusion voltage. The switch pulse had duration Tp = 0.1 /zs and the hold pulse had 
duration Thoid = 0.9 /j,s. The amplitude of the switch pulse ip was adjusted for each curve so that P(Ip) — 0.5. In 
this simulation it was first determined that a pulse consisting of the hold level itself (i.e. when ip = ihoid) did not 
induce any switching (P < 0.001). We can see from the solid curve of Fig. 0] that the hold level ihoU — 0.34 induces 
late switching events, with the actual switching taking place during the hold time, resulting in a measurement time 
that exceeds Tp. Lowering the hold level to itoid — 0.3 (dashed curve of Fig. 0} causes the switches to occur earlier 
in time, with a slight reduction of Tmeas- If the hold level is lowered to ihoid = 0.25 (dash-dotted curve of Fig.0J we 
find that the distribution P{t) becomes sharply peaked, with Tmeas ~ Tp. 

The simulations clearly show that Tmeas increases with ihoid- We find Tmeas = 0.10,0.42,0.44 /zs for ihoid = 
0.25,0.30,0.34 respectively. These values of Tmeas are approximately the same for different switching probabilities 



in the region 



± Ai/2. For ihoid > 0.34, Tmeas will increase, since the hold pulse alone will begin to induce 



switching events. For ihoid < 0.25, Tmeas will continue to be close to Tp, although for very low ihoid it will become 
hard to distinguish between switch and no-switching events. We have also run this simulation with lower strength of 
the noise term, r„ — 0.11, and we observe Tmeas ~ 0.1 — 0.15 /is for 0.40 < ihoid < 0.63. As expected, a weaker noise 
term will give a more rapid measurement time over a wider range of hold levels. 

From this simulation we can conclude the following: 1) Even in the presence of significant noise (r„ — 0.47) it is 
possible to find a switch pulse amplitude and hold level which give a rapid dynamical switching that occurs during 
the switch pulse without a significant number of late switches. 2) Simply adjusting the hold level for zero switching 
when ip = ihoid is not good enough to achieve optimal measurement time. Because one can not observe the escape 
process in the experiment, simulations of this kind are necessary to verify that the hold level is sufficiently low. 

Several switching curves P{I) of the type shown schematically in Fig.|2fb) were simulated for various values of the 
temperature and switching pulse parameters, as they were adjusted in the actual experiment. In Fig. |Slwe show a 
comparison between the simulated and measured resolution AI/Is^, as a function of the switch pulse duration Tp. 
Here, Igw is the value of the switch pulse amplitude which gives 50% switching probability, P{Ip) = 0.5. It is difficult 
to compare the simulated and experimental value of Isw, because in the experiment we measure the voltage pulse 
applied to the bias resistor. Dispersion in the twisted pairs, and parasitic capacitance of the sample mount, make it 
difficult to know the actual amplitude of the current pulse applied to the junction. 

The simulated values shown in Fig. [S] are calculated for two different noise tempertures. The value r„ — 0.005 
(or T„ = 30 mK for the particular junction parameters considered), corresponds to the measured temperature of 
the sample mount during the experiment. For this value of the noise temperature we see that the experimental and 
simulated values agree for longer Tp, while for shorter pulses the simulation gives much higher resolution than the 
experiment. We have adjusted the noise temperature, and we find reasonable agreement over the range of Tp studied 
for r„ = 0.47 {T — 3 K). This temperature seems excessively high, but note that in the simulation we have put Rt 
and i?2 at the same temperature. It would be more realistic to take Ri, at a higher, and i?2 at a lower temperature. In 
turn, this would increase inb and decrease j„2 in Eq. ^, causing a partial cancellation effect, giving an intermediate 
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FIG. 4: This shows P{t) during one pulse, for dt = 0.01 ^s, P = 0.50 and Vtrig = 35 ^V, with r„ = 0.47 and ihoid = 0.34 (solid 
line), ihoid = 0.30 (dashed line) and ihoid = 0.25 (dash-dotted line). 




FIG. 5: This shows AI /Ism as a function of Tp. The dashed curves are simulated data with r„ — 0.005 and 0.47, and the solid 
is measured data with F = 0.005 (T = 0.030 K). 



effective temperature. The exact consequences remains though to be investigated. 

Both the experimental and theoretical curves of Fig. \E\ show decreasing resolution (increasing A/) for shorter 
measurement time (smaller Tp). Partly this is due to that isw increases with decreasing Tp for given r„, but also the 
statistical fluctuation of the mean noise amplitude should be bigger for a short pulse than for a long pulse. Thus, as 
expected, there is a trade off between measurement speed and measurement resolution. 



V. SUMMARY 

In summary, we have performed computer simulations of the dynamical switching process of a Josephson junction 
subject to noise and to frequency dependent damping. The model for damping is appropriate for measurements of small 
capacitance Josephson junctions, where the phase dynamics is overdamped at high frequencies, and underdamped at 
low frequencies. We use a special bias current pulse, consisting of a short switch pulse followed by a longer hold level, 
appropriate for binary detection of the switching process. We characterize the switching process in terms of speed, 
Tmeas, and resolution, AI/Isw The simulations show that it is possible to achieve a minimum switching time given 
by the duration of the fast switching pulse, even in the presence of strong noise. To achieve this fast measurement, the 
hold level must be set appropriately low. As expected, there is a trade off between measurement speed and resolution. 
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